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Abstract 

Adiabatic quantum computation has recently attracted attention in the physics and com- 
puter science communities, but its computational power was unknown. We describe an effi- 
cient adiabatic simulation of any given quantum algorithm, which implies that the adiabatic 
computation model and the conventional quantum computation model are polynomially equiv- 
alent. Our result can be extended to the physically realistic setting of particles arranged on a 
two-dimensional grid with nearest neighbor interactions. The equivalence between the models 
provides a new vantage point from which to tackle the central issues in quantum computation, 
namely designing new quantum algorithms and constructing fault tolerant quantum computers. 
In particular, by translating the main open questions in the area of quantum algorithms to the 
language of spectral gaps of sparse matrices, the result makes these questions accessible to a 
wider scientific audience, acquainted with mathematical physics, expander theory and rapidly 
mixing Markov chains. 

1 Introduction 



The study of adiabatic quantum computation was initiated several years ago by Farhi, Goldstone, 
Gutmann and Sipser [T3j , who suggested a novel quantum algorithm for solving classical optimiza- 
tion problems such as Satisfiability (Sat). Their algorithm is based on a celebrated theorem 
in quantum mechanics known as the adiabatic theorem |19U25| . Although the exact worst-case 
behavior of this algorithm is not known, several simulations (see, e.g., ^3]) ° n random instances 
of up to 20 quantum bits led to various optimistic speculations. The bad news is that there is now 
mounting evidence pi llUll^] that the algorithm of |14| takes exponential time in the worst-case 
for NP-complete problems. Nevertheless, adiabatic computation was since shown to be promising 
in other less ambitious directions: it possesses several interesting algorithmic capabilities, as we 
will soon review, and in addition, it exhibits inherent robustness against certain types of quan- 
tum errors jH]. We note that a small scale adiabatic algorithm has already been implemented 
experimentally, using a Nuclear Magnetic Resonance (NMR) system |36j . 
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We briefly describe the model of adiabatic computation (a more precise description appears 
in Section I2.2|) . A computation in this model is specified by two Hamiltonians named -ffimt and 
-f^finai ( a Hamiltonian is simply a Hermitian matrix). The eigenvector with smallest eigenvalue (also 
known as the ground state) of H{ n [ t is required to be an easy to prepare state, such as a tensor 
product state. The output of the adiabatic computation is the ground state of the final Hamiltonian 
.fffinai. Hence, we choose an Hfi na i whose ground state represents the solution to our problem. We 
require the Hamiltonians to be local, i.e., we require them to only involve interactions between a 
constant number of particles (this can be seen as the equivalent of allowing gates operating on 
a constant number of qubits in the standard model). This, in particular, makes sure that the 
Hamiltonians have a short classical description, by simply listing the matrix entries of each local 
term. The running time of the adiabatic computation is determined by the minimal spectral gap 1 
of all the Hamiltonians on the straight line connecting H m \ t and H^ na \: H(s) = (1 — s)Hi n { t + si7g na i 
for s £ [0, 1]. More precisely, the adiabatic computation is polynomial time if this minimal spectral 
gap is at least inverse polynomial. 

The motivation for the above definition is physical. The Hamiltonian operator corresponds to 
the energy of the quantum system, and for it to be physically realistic and implementable it must 
be local. Its ground state is the state of lowest energy. We can set up a quantum system in the 
ground state of H m \ t (which is supposed to be easy to generate) and apply the Hamiltonian -ffmit 
to the system. We then slowly modify the Hamiltonian along the straight line from -ffjmt towards 
final- It follows from the adiabatic theorem that if this transformation is performed slowly enough 
(how slow is determined by the minimal spectral gap), the final state of the system will be in the 
ground state of i?fl na i, as required. 

What is the computational power of this model? In order to refer to the adiabatic model as a 
computational model that computes classical functions (rather than quantum states), we consider 
the result of the adiabatic computation to be the outcome of a measurement of one or more of 
the qubits, performed on the final ground state. It is known that adiabatic computation can be 
efficiently simulated by standard quantum computers [HIE]- Hence, its computational power is 
not greater than that of standard quantum computers. Several positive results are also known. 
In |9ll30| it was shown that Grover's quadratic speed-up for an unsorted search can be realized 
as an adiabatic computation. Moreover, |HII29ll32j showed that adiabatic computation can 'tunnel' 
through wide energy barriers and thus outperform simulated annealing, a classical counterpart of the 
adiabatic model. However, whether adiabatic computation can achieve the full power of quantum 
computation was not known. In fact, even the question of whether adiabatic computation can 
simulate general classical computations efficiently was unknown. The focus of this paper is the 
exact characterization of the computational power of adiabatic computation. 

Before we describe our results, let us clarify one subtle point. Most of the previous work on the 
subject focused on a restricted class of adiabatic algorithms that can be referred to as adiabatic 
optimization algorithms. In these algorithms, -fffinal is chosen to be a diagonal matrix, corresponding 
to a combinatorial optimization problem. In particular, this implies that the ground state of H^ a& \ 
(which is the output of the computation) is a classical state, i.e., a state in the computational 
basis. In this paper, however, we associate the term adiabatic computation with the more general 
class of adiabatic algorithms, where the only restriction on -£/fi na i is that it is a local Hamiltonian. 
We do this because, from a physical point of view, there is no reason to force the physical process 
described above to have a diagonal H g na i, when all other Hamiltonians are not restricted this way. 
Thus, our definition of adiabatic computation seems to be the natural one to use. It is this natural 
definition that allows us to prove our results. 

lr The spectral gap is the difference between the lowest and second lowest eigenvalue. 
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1.1 Results - Computational Complexity of the Adiabatic Model 

Our main result clarifies the question of the computational power of adiabatic algorithms. We 
show: 

Theorem 1.1 The model of adiabatic computation is polynomially equivalent to the standard model 
of quantum computation. 

As mentioned above, one direction of the equivalence is already known [13115] . Our contribution 
is to show that standard quantum computation can be efficiently simulated by adiabatic compu- 
tation. We do this by using adiabatic computation with 3-local Hamiltonians. We note that [I] 
made a preliminary step in the direction of Theorem 11.11 but the model that they considered was 
quite different. 2 

One corollary of our main theorem is the following. We can consider the model of adiabatic 
computation with a more general set of Hamiltonians known as explicit sparse Hamiltonians. These 
are Hermitian matrices that have at most polynomially many nonzero elements in each row and 
column, and, moreover, there is an efficient Turing machine that can generate a list of all nonzero 
entries in a given row or column. Clearly, local Hamiltonians are a special case of explicit sparse 
Hamiltonians. It was shown in [1] that adiabatic computation with explicit sparse Hamiltonians 
can still be simulated by standard quantum computation (this extends the result of j5JE] i n a 
non-trivial way). Hence, we obtain the following result. 

Corollary 1.2 The model of adiabatic computation with explicit sparse Hamiltonians is polynomi- 
ally equivalent to the standard model of quantum computation. 

Explicit sparse matrices are pervasive in computer science and combinatorics, and hence this corol- 
lary might be more useful than Theorem 11.11 in the context of the design of quantum algorithms 
and the study of quantum complexity. 

To summarize, our results show that questions about quantum computation can be equivalently 
considered in the model of adiabatic computation, a model that is quite different from the more 
common circuit-based models. There are two reasons why it seems worthwhile to try to design 
quantum algorithms in the adiabatic framework. First, there are several known powerful techniques 
to analyze spectral gaps of matrices, including expander theory ^3] and rapidly mixing Markov 
chains 24,34 . Indeed, probability theory is often used in mathematical physics to analyze spectral 
gaps of Hamiltonians (see, e.g., |3S])) and our proofs also make extensive use of Markov chain tools. 
Second, it is known that many interesting algorithmic problems in quantum computation can be 
cast as quantum state generation problems [I]. The problem of generating special quantum states 
seems more natural in the adiabatic model than in the standard model. 

1.2 Results - Towards Experimental Implications 

Theorem 1 1 . 1 1 uses 3- local Hamiltonians that act on particles that may be arbitrarily far apart. From 
a practical point of view, it is often difficult to create controlled interactions between particles 
located far-away from each other. Moreover, three-particle Hamiltonians are technologically very 
difficult to realize. If one wants to physically realize the adiabatic algorithms, it would be much 
better to have only two-local interactions between nearest neighbor particles. To this end we prove 

2 Namely, |1] showed that adiabatic computation using simulatable Hamiltonians is as powerful as standard quan- 
tum computation. Simulatable Hamiltonians are Hamiltonians that can be simulated efficiently by a quantum circuit. 
They are very different from local Hamiltonians, and they cannot even be written explicitly. Instead, such Hamilto- 
nians are specified using products of local unitary matrices. 
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the following theorem. This, we believe, brings the adiabatic computation model one step closer to 
physical realization. 

Theorem 1.3 Any quantum computation can be efficiently simulated by an adiabatic computation 
with two-local nearest neighbor Hamiltonians operating on six-state particles set on a two dimen- 
sional grid. 

The need for six-state particles arises from our construction. It is an open question whether this 
can be improved. 

Theorems 11.11 and 11.31 open up the possibility of physically realizing universal quantum com- 
putation using adiabatically evolving quantum systems. As mentioned before, there are possible 
advantages to this approach: adiabatic quantum computation is resilient to certain types of noise |Hj. 
An important component of this resilience is the existence of a spectral gap in the Hamiltonian. 
It is well known in physics that such a gap plays an important role in the context of protecting 
quantum systems from noise. However, it remains to be further studied, both experimentally and 
theoretically, what the right model for noisy adiabatic computation is, and whether fault tolerant 
adiabatic computation can be achieved. We refer the reader to further discussion in Subsection II .61 

1.3 Proof of Theorem II. H Overview 

Given an arbitrary quantum circuit |26j . our goal is to design an adiabatic computation whose 
output is the same as that of the quantum circuit. Some similarities between the models are 
obvious: one model involves unitary gates on a constant number of qubits, while the other involves 
local Hamiltonians. However, after some thought, one eventually arrives at the following difficulty 
The output state of the adiabatic computation is the ground state of H g na i. The output state of 
the quantum circuit is its final state, which is unknown to us. How can we specify -fffinai without 
knowing the output state of the quantum circuit? Notice that this state can be some complicated 
quantum superposition. One might wonder why our task is not trivial, since this state does have 
an efficient local classical description, namely the quantum circuit. However, local quantum gates, 
which operate in sequence to generate a non-local overall action, are very different from local 
Hamiltonians, which correspond to simultaneous local constraints. To explain the solution, we first 
set some notations. 

Without loss of generality we assume that the input to the quantum circuit consists of n qubits 
all initialized to |0)'s. 3 Then, a sequence of L unitary gates, U\, . . . ,Ul, each operating on one or 
two qubits, is applied to the state. The system's state after the fth gate is \a(£)). The output 
of the quantum circuit is in general a complicated quantum state \a{L)) of n qubits, which is 
then measured in the standard basis. We now want to associate with it a corresponding adiabatic 
computation. 

A first natural attempt would be to define H g na i as a local Hamiltonian with \a{L)) as its ground 
state. However, this attempt encounters the difficulty mentioned above: not knowing \a(L)), it 
seems impossible to explicitly specify //final- The key to resolve this difficulty is the observation that 
the ground state of -//final need not necessarily be the state \a(L)). It is sufficient (under some mild 
restrictions) that the ground state has a non-negligible inner product with \a(L)). This gives us 
significant flexibility in designing H^ na \. Our idea is to base our solution on a seemingly unrelated 
ingenious result of Kitaev [22], in which he provides the first quantum NP-complete problem, 
namely, LOCAL Hamiltonians. This result can be viewed as the quantum analogue of the Cook- 
Levin theorem [2S], which states that 3-Sat is NP-complete. For his proof, Kitaev defined a local 

3 Otherwise, the first n gates can be used to flip the qubits to the desired input. 
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Hamiltonian that checks the time propagation of a quantum circuit. Kitaev's local Hamiltonian 
has as its ground state the entire history of the quantum computation, in superposition: 



The right {L qubits) register is a clock that counts the steps by adding Is from left to right. The 
superscript c denotes clock qubits. We note that this state has a non-negligible projection on our 
desired state \a{L)). Hence, instead of designing a Hamiltonian that has the final unknown state 
of the circuit as its ground state, a task that seems impossible, we can define -fffinai to be Kitaev's 
local Hamiltonian. Why is it possible to define a local Hamiltonian whose ground state is |ry), 
whereas the same task seems impossible with |a(L))? The idea is that the unary representation of 
the clock enables a local verification of correct propagation of the computation from one step to 
the next, which cannot be done without the intermediate computational steps. 

We thus choose Kitaev's Hamiltonian 22. to be our -fffi na i. This Hamiltonian involves five body 
interactions (three clock particles and two computation particles). For the initial Hamiltonian 
Hi n \t we require that it has |a(0)) (g) |0 L ) C , the first term in the history state, as its unique ground 
state. It is easy to define such a local Hamiltonian, because |a(0)) (8> |0 L ) C is a tensor product 
state. Crucially, H iri \ t and -fffmai can be constructed efficiently from the given quantum circuit; no 
knowledge of \a{L)) is required for the construction. 

A technical problem lies in showing that the spectral gap of the intermediate Hamiltonian H (s) 
is lower bounded by some inverse polynomial (more specifically, we show it is larger than 1/L 2 ). To 
do this, we use a mapping of the Hamiltonian to a Markov chain corresponding to a random walk 
on the L + 1 time steps. We then apply the conductance bound from the theory of rapidly mixing 
Markov chains |34| to bound the spectral gap of this chain. We note that, in general, applying the 
conductance bound requires knowing the limiting distribution of the chain, which in our case is 
hard since it corresponds to knowing the coefficients of the ground state for all the Hamiltonians 
H(s). We circumvent this problem by noticing that it is actually sufficient in our case to know very 
little about the limiting distribution of the Markov chain, namely that it is monotone (in a certain 
sense to be defined later). This allows us to apply the conductance bound, and deduce that the 
spectral gap is 0(1/L 2 ). From this it follows that the running time of the adiabatic computation is 
polynomial. Extracting the output of the quantum circuit from the history state efficiently is easy: 
Measure all the qubits of the clock and if the clock is in the state the computational qubits 
carry the result of the circuit. Otherwise, start from scratch. 4 

The above scheme gives a proof of Theorem 11.11 that uses 5-local Hamiltonians, and runs in 
time roughly 0(L 5 ). The improvement to 3-locality is based on a simple idea (used in [20 to prove 
that the 3-local Hamiltonian problem is quantum NP-complete) but obtaining a lower bound on 
the spectral gap is significantly more involved technically. We postpone its explanation to the body 
of the paper. The running time we achieve in this case is roughly 0(L 14 ). 

1.4 Proof of Theorem II. 31 Overview 

The idea underlying the proof of Theorem 11.11 by itself does not suffice to prove Theorem 11.31 
The basic problem lies in arranging sufficient interaction between the computational and the clock 

4 This gives an overhead factor of L which can be avoided by adding O(^L) identity gates to the quantum circuit 
at the end, which has the effect that most of the history state is concentrated on the final state \a(L)). See 
S ubsect ion IM . 1 , M I for more details. 




(1) 
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particles, since if the particles are set on a grid, each clock particle can only interact with four neigh- 
bors. We circumvent this problem as follows. Instead of having separate clock and computational 
particles, we now assign to each particle both clock and computational degrees of freedom (this is 
what makes our particles six-state). We then construct a computation that propagates locally over 
the entire set of particles, snaking up and down each column of the lattice. The adiabatic evolution 
now ends up in the history state of this snake-like sequence of states. 

The lower bound on the spectral gap is obtained in an essentially identical way as in the 3-local 
Hamiltonian case. 

1.5 Related Work 

After the preliminary version of this paper appeared [2j, the results regarding Quantum-NP com- 
pleteness were tightened by [2L to show that the 2-local Hamiltonian problem is Quantum-NP 
complete. Following the ideas presented in the current paper, |2j used their result to show that 
Theorem II .11 holds when the Hamiltonians are 2-local. 

The idea to use an inverse polynomial spectral gap for fault tolerance is certainly not new. It is 
a crucial ingredient in topological (and later, geometrical) quantum computation |18|l23|l2*7] . Note, 
however, that in those models the spectral gap has no effect on the running time or on any other 
algorithmic aspects, and it is used only to separate the computational subspace from the "noisy" 
subspace. In contrast, the spectral gap in adiabatic computation is crucial from the algorithmic 
point of view, since it determines the time complexity of the computation. 

1.6 Open Questions 

This paper demonstrates that quantum computation can be studied and implemented entirely 
within the adiabatic computation model, without losing its computational power. This result 
raises many open questions in various directions. First, it would be interesting to determine if the 
parameters presented in this work can be improved. For example, it might be possible to shorten 
the running time of our adiabatic simulation. Decreasing the dimensionality of the particles used 
in Theorem II ..SI from six to two or three might be important for implementation applications. An 
interesting question is whether Theorem 11.31 can be achieved using a one dimensional instead of a 
two dimensional grid. 

Second, the possibility of fault tolerant adiabatic computation deserves to be studied both 
experimentally and theoretically. Since the publication of the preliminary version of the current 
paper [2|, several researchers have begun to study adiabatic computation in the presence of noise 
[Tl l31ll3*3*] . However, it is still unclear whether adiabatic evolution might be helpful for the physical 
implementation of quantum computers. 

Our results imply the equivalence between standard quantum computation and various other 
variants of adiabatic computation that have been considered in the literature and that are more 
general than our model. These include adiabatic computation with a general path between Hi n [ t 
and -ff final) rather than a straight line (see [3] and ^2] for a rigorous definition), and adiabatic 
computation with explicit sparse Hamiltonians [1] (see Corollary ll.2|) . A problem we leave open 
is to characterize the computational power of adiabatic optimization, studied in [9lllD|[T3]. In this 
model, the initial state is a tensor product of qubits in the state ^j(|0) + the final Hamiltonian 
is diagonal, and the evolution is carried out on a straight line. It is still possible that such a 
computation can be simulated efficiently by a classical Turing Machine. 

Finally, we hope that the adiabatic framework might lead to the discovery of new quantum 
algorithms. As shown in this paper, as well as in 0j, tools from probability theory, mathematical 
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physics and spectral gap analysis might turn out to be relevant and useful. In order to improve 
our understanding of the benefits of the adiabatic paradigm, it might be insightful to see adiabatic 
versions of known quantum algorithms, presented in a meaningful way. 

Organization: In Section [2] we describe the model of adiabatic computation and state some 
relevant facts about Markov chains. Section |3] shows how adiabatic systems with local Hamilto- 
nians allowing five- and later three-body interactions, can efficiently simulate standard quantum 
computations. Section |1] shows how to adapt the construction to a two-dimensional grid. 

2 Preliminaries 

2.1 Hamiltonians of n-Particle Systems 

For background on n-qubit systems, quantum circuits and Hamiltonians, see |26j . An n-particle sys- 
tem is described by a state in Hilbert space of dimension d n , the tensor product of n <i-dimensional 
Hilbert spaces. For simplicity, we restrict our discussion in this subsection to quantum systems 
composed of 2-dimensional particles, i.e., qubits; a similar discussion holds for higher dimensional 
particles (such as the 6-dimensional case we consider later). 

In the standard model of quantum computation, the state of n qubits evolves in discrete time 
steps by unitary operations. In fact, the underlying physical description of this evolution is con- 
tinuous, and is governed by Schrodinger's equation: — i-^\ip(t)) = H(t)\ip(t)). Here \ip(t)) is the 
state of the n qubits at time t, and H(t) is a Hermitian 2 n x 2 n matrix operating on the space 
of n qubits. This H{t) is the Hamiltonian operating on a system; it governs the dynamics of the 
system. Given that the state of the system at time t = is equal to IV'(O)), one can in principle 
solve Schrodinger's equation with this initial condition, to get \tjj(T)), the state of the system at a 
later time t = T. The fact that the Hamiltonian is Hermitian corresponds to the familiar fact that 
the discrete time evolution of the quantum state from time t± to a later time £2 is unitary. 

We sometimes refer to eigenvalues of Hamiltonians as energies. The ground energy of a Hamil- 
tonian is its lowest eigenvalue and the corresponding eigenvector (s) are called ground state(s). We 
define A(iJ), the spectral gap of a Hamiltonian H, to be the difference between the lowest eigen- 
value of H and its second lowest eigenvalue. (A(H) = if the lowest eigenvalue is degenerate, 
namely, has more than one eigenvector associated with it). We define the restriction of H to some 
subspace S, denoted H$, as IlgHTls where Hs is the orthogonal projection on S. 

A Hamiltonian on an n-particle system represents a certain physical operation that one can, in 
principle, apply to an n-particle system. However, it is clear that one cannot efficiently apply any 
arbitrary Hamiltonian (just describing it requires roughly 2 2n space). We say that a Hamiltonian 
H is k-local if H can be written as Y^A H A where A runs over all subsets of k particles, and H A 
operates trivially on all but the particles in A (i.e., it is a tensor product of a Hamiltonian on A 
with identity on the particles outside of A). Notice that for any constant k, a fc-local Hamiltonian 
on n-qubits can be described by 2 2k n k = poly(n) numbers. We say that H is local if H is /c-local 
for some constant k. 

In this paper we restrict our attention to fc-local Hamiltonians. This requirement corresponds to 
the fact that all known interactions in nature involve a constant number of particles. We attempt 
to make k as small as possible to make the Hamiltonian easier to implement. 



7 



2.2 The Model of Adiabatic Computation 



The cornerstone of the adiabatic model of computation is the celebrated adiabatic theorem |19U25j . 
Consider a time-dependent Hamiltonian H(s), s E [0,1], and a system initialized at time t = 
in the ground state of H(0) (here and in the following we assume that for all s E [0, 1], H(s) has 
a unique ground state). Let the system evolve according to the Hamiltonian H(t/T) from time 
t = to time T. We refer to such a process as an adiabatic evolution according to H for time 
T. The adiabatic theorem affirms that for large enough T the final state of the system is very 
close to the ground state of H(l). Just how large T should be for this to happen is determined by 
the spectral gap of the Hamiltonians H(s). Such an upper bound on T is given in the following 
theorem, adapted from [55] (whose proof in turn is based on see also [3] for a recent elementary 
proof of a slightly weaker version) . 

Theorem 2.1 (The Adiabatic Theorem (adapted from [29])) Let H\ n \ t and Hfi na \ be two Hamil- 
tonians acting on a quantum system and consider the time- dependent Hamiltonian H(s) ■= (1 — 
s)H{ n n + s.Hfi na i. Assume that for all s, H(s) has a unique ground state. Then for any fixed 5 > 0, 



then the final state of an adiabatic evolution according to H for time T (with an appropriate setting 
of global phase) is e-close in £2-norm to the ground state of H^ n . A \. The matrix norm is the spectral 
norm \\H\\ := max ffi ||iJu;[|/||u;||. 

One should think of 5 as being some fixed constant, say, 0.1. We cannot take 5 = because of the 
constant hidden in the f2 notation, which goes to infinity as 8 goes to 0. 

Let us now describe the model of adiabatic computation. In this paper we use the following 
definition of adiabatic computation that slightly generalizes that of Farhi et al. ^3]. The adiabatic 
'circuit' is determined by Hi n it and -fffi na i and the output of the computation is (close to) the ground 

State of -fffinal- 

Definition 2.2 A k-local adiabatic computation AC (n, d, -ffi n it> -Hflnab e ) is specified by two k-local 
Hamiltonians, H[ n a and -£ffi na i acting on n d- dimensional particles, such that both Hamiltonians 
have unique ground states. The ground state of H[ n n is a tensor product state. The output is a 
state that is e-close in li-norm to the ground state of H^ n& \. Let T be the smallest time such that 
the final state of an adiabatic evolution according to H{s) := (1 — s)H{ n n + signal for time T is 
e-close in £2-norm to the ground state of H^ na \. The running time of the adiabatic algorithm is 
defined to be T ■ max s [|i?(s)||. 

Observe that we have chosen our definition of running time to be T • max s [|iT(s)|| and not T. 
Notice that if the Hamiltonians are multiplied by some factor, this divides the bound of Equation 
121 and hence T, by the same factor. Hence, if the running time is defined to be T one would be able 
to achieve arbitrarily small running times, by multiplying the Hamiltonians by large factors. Our 
definition, on the other hand, is invariant under a multiplication of the Hamiltonian by an overall 
factor, and so takes into account the known physical trade-off between time and energy. 5 

5 This trade-off between time and the norm of the Hamiltonian (namely, the energy), is manifested in Schrodinger's 
equation whose solution does not change if time is divided by some factor and at the same time the Hamiltonian is 
multiplied by the same factor. 



T 




(2) 
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The right hand side of Equation [21 can be used to provide an upper bound on the running time 
of an adiabatic computation. Hence, in order to show that an adiabatic algorithm is efficient, it is 
enough to use Hamiltonians of at most poly(n) norm, and show that for all s £ [0, 1] the spectral 
gap A(H(s)) is at least inverse polynomial in n. 

We note that in certain cases, it is possible to obtain a stronger upper bound on the running 
time. Indeed, assume there exists a subspace S such that for all s S [0, 1], H(s) leaves S invariant, 
i.e., H(s)(S) C S. Equivalently, H(s) is block diagonal in S and its orthogonal space S- 1 . Consider 
Hs(s), the restriction of H(s) to S. Then, starting from a state inside S, an adiabatic evolution 
according to H is identical to an adiabatic evolution according to H$ (this follows from Schrddinger's 
equation). Hence, we can potentially obtain a stronger upper bound by replacing A(H(s)) with 
A(Hg(s)) in Equation |2J This stronger upper bound will be used in our first adiabatic algorithm. 

Finally, let us mention that one can define more general models of adiabatic computation. For 
example, one might consider non-local Hamiltonians (see [I]). Another possible extension is to 
consider more general paths between fli n it and -fffi na i (see, e.g., pSIBllI])- Obviously, our main 
results, such as Theorem ll.il hold also for these more general models. 



2.3 Markov Chains and Hermitian Matrices 

Under certain conditions, there exists a standard mapping of Hamiltonians to Markov chains (for 
background on Markov chains, see |24j ) . The following fact is useful to show that this mapping 
applies in the case we analyze. 

Fact 2.3 (Adapted from Perron's Theorem, Theorem 8.2.11 in |17j ) LetG be a Hermitian 
matrix with real non-negative entries. If there exists a finite k such that all entries of G k are posi- 
tive, then G 's largest eigenvalue is positive, and all other eigenvalues are strictly smaller in absolute 
value. Moreover, the corresponding eigenvector is unique, and all its entries are positive. 

We define the mapping for G, a Hermitian matrix operating on an L + 1 dimensional Hilbert 
space. Suppose that all the entries of G are real and non-negative, that its eigenvector (ao, ...,«£,) 
with largest eigenvalue \i satisfies on > for all < i < L and that fj, > 0. Define P by: 

Pij := ~Gij. (3) 

The matrix P is well defined, and is stochastic because all its entries are non-negative and each of 
its rows sums up to one. It is easy to verify the following fact: 

Fact 2.4 The vector (vq, ■ ■ ■ , v£) is an eigenvector ofG with eigenvalue 5 if and only if(atQVo, ■ ■ ■ , q?l^l) 
is a left eigenvector of P with eigenvalue S/fi. 

We will consider G of the form G = I — H for some Hamiltonian H. The above fact implies that 
if (ao, . . . , ai) is the ground state of H with eigenvalue A then (oq, . . . , a 2 L ) is a left eigenvector 
of P with maximal eigenvalue 1. By normalizing, we obtain that it ■= (a^/Z, . . . , a\/Z) is the 
limiting distribution of P, where Z = ^af. Moreover, the gap between P's largest and second 
largest eigenvalues is equal to A(H)/(1 — A). 

2.4 Spectral Gaps of Markov Chains 

Given a stochastic matrix P with limiting distribution it, and a subset B C {0, . . . ,L}, the flow 
from B is given by F(B) := YlieBj&B ^i^ij- Define the 7r-weight of B as n(B) := J2ieB n i- ^ ne 
conductance of P is defined by <p(P) '■= min^ F(B)/tt(B), where we minimize over all non-empty 
subsets B C {0, . . . , L} with n(B) < \. 
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Theorem 2.5 (The conductance bound |34j)) The eigenvalue gap of P is at least ^ip(P) 2 . 

3 Equivalence of Adiabatic and Quantum Computation 

Here we prove Theorem II .!( by showing how to simulate a quantum circuit with L two-qubit gates 
on n qubits by an adiabatic computation onn + L qubits (the other direction was shown in |14ll§] ) . 
We first allow five qubit interactions. We later show how to reduce it to three, using techniques 
that will also be used in Section |IJ 

3.1 Five-local Hamiltonian 

Theorem 3.1 Given a quantum circuit on n qubits with L two-qubit gates implementing a uni- 
tary U, and e > 0, there exists a 5-local adiabatic computation AC(n + L, 2, H m it, -fffmab e) whose 
running time is poly(L, -) and whose output (after tracing out some ancilla qubits) is e-close (in 
trace distance) to U\0 n ). Moreover, H[ n i t and H^ n& \ can be computed by a polynomial time Turing 
machine. 

The running time we obtain here is 0(e~( 5 + 3<5 )L 5+2<5 ) for any fixed 5 > 0. 
3.1.1 The Hamiltonian 

For our construction we use the Hamiltonian defined in |22j . Denote \je) := \a(£)) ® \l i L ~ e ) c , 
where \a{tj) denotes the state of the circuit after the ith. gate and the superscript c denotes the 
clock qubits. We would like to define a local Hamiltonian H in a with ground state I70) = \0 n ) ®\0 L } C , 
and a local Hamiltonian H^ n& \ with ground state \rj) = Yld=o \lt) as m Equation ^ To do 

this, we write flmit and H^ n& \ as a sum of terms: 

-f^init • -^clockinit ~t~ -^input -^clock 

-f^final : = - Hi + Hinput + -^clock- 
1=1 

The terms in H^ n& \ (and likewise in H^ n \ t ) are defined such that the only state whose energy (i.e., 
eigenvalue) is is the desired ground state. This is done by assigning an energy penalty to any 
state that does not satisfy the required properties of the ground state. The different terms, which 
correspond to different properties of the ground states, are described in the following paragraphs. 
The adiabatic evolution then follows the time-dependent Hamiltonian 

H(s) = {l-s)H init + sH^ h (4) 

Notice that as s goes from to 1, i? c iockinit is slowly replaced by \ J2e=i He while H mput and -ff c iock 
are held constant. 

We now describe each of the terms. First, -ffdock checks that the clock's state is of the form 
\l e L ~^) c for some < i < L. This is achieved by assigning an energy penalty to any basis state 
on the clock qubits that contains the sequence 01, 

L-l 

#clock := El 01 ><°%+1> 

1=1 
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where the subscript indicates which clock qubits the projection operates on. Note that illegal clock 
states are eigenstates of -ff c iock with eigenvalue at least 1; legal clock states have eigenvalue 0. 
Next, H mput checks that if the clock is |0 L ) C , the computation qubits must be in the state |0 ra ), 



tfinput :=X)|l><l|i® |0)(0|5. 



i=l 

We complete the description of -ff; n it with -ff c iockinit whose goal is to check that the clock's state 
is |0 L ) C , 

-^clockinit := 

Claim 3.2 The state \jq) is a ground state of H m \t with eigenvalue 0. 6 

Proof: It is easy to verify that -ffmit|7o) =0. As a sum of projectors, -ffmit is positive semidefmite 
and hence I70) is a ground state of H m \ t . ■ 

We now proceed to the first term in iJfl na i- The Hamiltonian Hi checks that the propagation 
from step £ — 1 to I is correct, i.e., that it corresponds to the application of the gate Ug. For 
1 < £ < L, it is defined as 

H t := I®|100)(100||_ w+1 -^(»|110)(100|t w+ i 

-U} ® |100)(110|t w+ i +/<E» |110)(110|ti )W - (5) 



Intuitively, the three-qubit terms above move the state of the clock one step forward, one step 
backward, or leave it unchanged. The accompanying matrices Ui, lit describe the associated time 
evolution. For the boundary cases £ = 1,L, we omit one clock qubit from these terms and define 



E x := 1(8) |00)(00|i, 2 - Ux ® |10)(00|i, 2 - u\ ® |00)(10|i, 2 + 1 ® |10)(10|i, 2 

H L := I®\1Q){1Q\l-x,l-U l ® \11){1Q\l-i,l-U[® |10)<11U_i,l + J® (6) 

Claim 3.3 The history state \rj) is a ground state of -fffinai with eigenvalue 0. 

Proof: It is easy to verify that H^ na \\rj) = 0. It remains to notice that for all 1 < I < L, Hi is 
positive semidefmite and hence so is Hfl na i. ■ 

Remark: Strictly speaking, Theorem 13 . 1 1 holds even if we remove the terms -ff c iock and -Hmput from 
both -ffmit and H fi na j. We include them mainly for consistency with the rest of the paper. 



3.1.2 Spectral Gap in a Subspace 

Let So be the L + 1-dimensional subspace spanned by I70), . . . , \ It is easy to verify the following 
claim. 

Claim 3.4 The subspace Sq is invariant under H(s), i.e., H(s){Sq) C Sq 

In this subsection, we show that the spectral gap of Hg Q (s), the restriction of H to So, is inverse 
polynomial in L. As mentioned in Subsection XL .21 this, together with Claim|23J is enough to obtain 
a bound on the running time of the adiabatic algorithm. 

6 The state I70) is in fact the unique ground state of -Hinit as will become apparent from the proof of the spectral 
gap. A similar statement holds for Claim HP1 
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Lemma 3.5 The spectral gap of the restriction of H(s) to Sq satisfies A(Hs (s)) = Q(L 2 ) for all 

se[0,l}. 



Proof: Let us write the Hamiltonians i?s 0j init and fls 0) final in the basis I70), . . . , |7l) of Sq. Both 
-f^ciock and -ffinput are on Sq and can thus be ignored. We have the following [L + 1) x (L + 1) 
matrices: 



/ ... \ 
1 ... 



\0 ... 1/ 



(7) 



Hs , final 



5I70X70I - ||to><Ti| - iI7l)(7l-i| + ^7l)(7z 



L-l 



/ 3 

l 

2 



\ 



V 




1 



(8) 



2 2 / 



We now lower bound A(Hs (s)). We consider two cases: 

The case s < 1/3: Here, Hs Q {s) is sufficiently close to Hs 0) init (whose spectral gap is 1) so we can 
apply the following standard lemma (see, e.g., [7], Page 244). 

Lemma 3.6 (Gerschgorin's Circle Theorem) Let A be any matrix with entries a«. Consider 
the discs in the complex plane given by 



D, 



\z - au\ < 



1 < i < n. 



Then the eigenvalues of A are contained in U-Dj and any connected component of UDj contains as 
many eigenvalues of A as the number of discs that form this component. 

For s < 1/3, Hso(s)i t i < 1/6 and Ylj^i H s ( s )i,j < V 6 - Moreover, for any i / 1, H So (s) iti > 
5/6 and Ylj^i Hso( s )i,j < 1/6- By the above lemma, we obtain that there is one eigenvalue smaller 
than 1/3 while all other eigenvalues are larger than 2/3. Hence, the spectral gap is at least 1/3. 

The case s > 1/3: We note that Hs final is t ne Laplacian of the simple random walk [25 of a 
particle on a line of length L + 1. A standard result in Markov chain theory implies A(H$ 0t final) — 
0,(1/1?) [2U- For s > 1/3, H So (s) is sufficiently close to H So , ,final to apply Markov chain techniques, 
as we show next. 

Let (ao, • • • , ul)^ be the ground state of Hs (s) with eigenvalue A. Define the Hermitian matrix 
G(s) = I — Hs (s). It is easy to see that G(s) satisfies the conditions of Fact 12.31 for all s > 0. 
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We obtain that the largest eigenvalue \i = 1 — A of G(s) is positive and non-degenerate and the 
corresponding eigenvector («o, • • • , has positive entries. We can now map the matrix G(s) to 
a stochastic matrix P(s) as described in Subsection 12.31 The transition matrix P(s) describes a 
random walk on the line of L + 1 sites (Fig. Fact 12.41 implies that the limiting distribution of 
P(s) is given by ir = (oq/Z, . . . , a 2 L /Z) where Z = af. 

Pk,k-1 Pk,k+1 

o — o — o — o— o— o — o — o 

1 ... k - 1 k k + l ... L 

Figure 1: The random walk of P(s) 

We bound the spectral gap of P(s) using the conductance bound (see Subsection 12. 4|) . To do 
this we need to know that ir is monotone. We first show: 

Claim 3.7 For all < s < 1, the ground state of Hs (s) is monotone, namely ao > a± > . . . > 
OL L > 0. 

Proof: The case s = is obvious, so assume s > 0. We first claim that the ground state 
(ao, • • • , a.£y of Hs (s) = I — G(s) can be written as the limit 

- lim(G( S )/^(l,...,l)t 

Co I— >oo 

for some constant cq > 0. To see this, let \vq), . . . , \vl) be an orthonormal set of eigenvectors of 
G(s), with corresponding eigenvalues (J® > fii > . ■ . > (jll- By Fact 12 .31 the largest eigenvalue 
corresponds to a unique eigenvector, and hence we have \vq) = (ao, . . . , chl), and fiQ = fi. 

The set of eigenvectors \vi) forms an orthonormal basis, and we can write (1, . . . , l)t in terms 
of this basis: (1, . . . , l)t = ^ a\vi). Now, we have that l) f = Ei CiifYlvi)- By 

Fact 12 .31 we have < /i for all « 7^ 0, and \i > 0. We thus have that lim^ 00 (G(s)//i) £ (l, . . . , ly = 
co\v ). 

It is easy to check that G(s) preserves monotonicity, namely, if G(s) is applied to a monotone 
vector, the result is a monotone vector. Hence, when G(s)/n is applied to the monotone vector 
(1, ...,1)t, the result is a monotone vector. Thus, co|^o) is monotone. Finally, we observe that 
Co > 0. This is because Co is the inner product between the all 1 vector, and \v$) whose entries are 
all positive by Fact 12.31 This implies that \vq) is also monotone, as desired. ■ 

It follows that 7r is also monotone. We use this and simple combinatorial arguments to prove 
the following claim. 

Claim 3.8 For all 1/3 < s < 1, tp(P(s)) > ^. 

Proof: We show that for any nonempty B C {0, . . . , L}, F(B)/ir(B) > We consider two cases. 
First, assume that G B. Let k be the smallest such that k € B but k + l ^ B. Then, 



F{B) > ir k P(s) k k+1 = ir k ■ v^fc+i _ V^^k+i > lp?±LG(s) k k+1 

HyfKk 1 - A 1 - A 

where the last inequality follows from the monotonicity of ir. Using the definition of G and the 
assumption that s > 1/3 we get that G(s) kk +% > 1/6. We also have < 1 — A < 1, where the 
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second inequality follows from the fact that Hs (s) is positive semidefinite, and the first follows 
from fj, > which we previously deduced from Fact 12.31 Hence, 



*W > (9) 

By 7r(B) < 1/2, we have n{{k + 1, ... ,L}) > 1/2. Together with vr({fc + 1, . . . ,L}) < L-n k+1 we 
obtain vr fc+ i > 1/(2L). This yields the desired bound F(B)/n(B) > 1/(6L). 

Now assume that ^ B and let be the smallest such that k ^ B and A; + 1 G B. It is easy to 
see that irkP( s )k,k+i = ^fc+i-P^fc+i.fc- Hence, using the same argument as before we can see that 
Equation |§] holds in this case too. Since BC{H1,..., L}, we have ir({k + 1, . . . , L}) > n(B). 
Hence, ir k+1 > tt(B)/L. Again, this yields the bound F(B)/ir(B) > 1/(6L). ■ 

By Theorem 12, 5| we have that the spectral gap of P(s) is larger than 1/(2 • (6) 2 • L 2 ). By Sub- 
section [TJl we have that A(Hs ) > • (6) 2 L 2 ). Finally, notice that /j, = 1 — A > |> because 
A<(7o|i?5 (s)l7o) = |<i ■ 



3.1.3 Running Time 

We now complete the proof of Theorem 13.11 Note that we have already proved something which is 
very close to Theorem 13. II 

Claim 3.9 Given a quantum circuit on n qubits with L gates, the adiabatic algorithm with H{ n [ t 
and -fffmai as defined in the previous section, with T = 0(e- s L A+2S ) for some fixed 5 > 0, outputs 
a final state that is within ^-distance e of the history state of the circuit, \rj). The running time of 
the algorithm is 0(T ■ L). 

Proof: Claim E3 shows that So is invariant under H. Hence, as mentioned in Subsection 12.21 an 
adiabatic evolution according to H is identical to an adiabatic evolution according to H$ . Using 
Lemma |3.5I and Theorem 12.11 (with 1 1 -ffinlt — ^flnalll = 0(1)), we obtain that for T as above the 
final state (with global phase adjusted appropriately) is indeed e-close in ^2- norm to \n). By our 
definition, the running time of the adiabatic algorithm is 0(T ■ L) since ||i?(s)|| < (1 — s)||i^i n it|| + 
s || -ff final || = 0(L + n) = 0(L). The last equality follows from n = O(L), because each qubit is 
assumed to participate in the computation (otherwise we can omit it). ■ 

In fact, one might be satisfied with this claim, which enables generating adiabatically a state 
which is very close to \rj), instead of our desired \a(L)}. To see why this might be sufficient to 
simulate quantum circuits, suppose for a moment that e is 0, and the final state is exactly \r/). 
As mentioned in the introduction, we can now measure the clock qubits of the history state, and 
with probability 1/L the outcome is I = L, which means that the state of the first register is the 
desired state \a(L)). If the measurement yields another value, we repeat the adiabatic algorithm 
from scratch. To get I = L with sufficiently high probability, we repeat the process O(L) times, 
which introduces an overhead factor of L. The above discussion is also true with e > 0, as long as 
it is much smaller than 1/L, the weight of \a{L)) in \rj). 

However, this is not sufficient to complete the proof of Theorem 13.11 Indeed, the theorem as 
stated follows our definition of the model of adiabatic computation, which allows to perform one 
adiabatic evolution and then measure (and possibly trace out some qubits). Classical postprocessing 
such as conditioning on I being equal to L, and repeating the computation if it is not, are not 
allowed. Hence, we need to adiabatically generate a state that is close to \a(L)). 
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This technical issue can be resolved with the following simple trick, which at the same time 
allows us to avoid the overhead factor of L introduced before. We simply add another 0(-L) 
identity gates to the original quantum circuit at the end of its computation and then apply the 
adiabatic simulation to this modified circuit. This modification increases the weight of \a{L)) in 
the history state. The following lemma makes this precise. 

Lemma 3.10 Assume we can transform any given quantum circuit with L two-qubit gates on n 
qubits into a k-local adiabatic computation on n+L d- dimensional particles whose output is e close 
in (,2-norm to the history state of the quantum circuit and whose running time is f(L, e) for some 
function f . Then, we can transform any given quantum circuit with L two-qubit gates on n qubits 
into a k-local adiabatic computation on n+2L/e d-dimensional particles whose output (after tracing 
out some ancilla qubits) is e close in trace distance to the final state of the circuit and whose running 
time is /(2L/e,e/2). 

Proof: Given a quantum circuit on n qubits with L gates, consider the circuit obtained by ap- 
pending to it (~ — 1)L identity gates. Let L' = 2L/e be the number of gates in the modified circuit 
and let \rj) denote its history state. By our assumption, we can transform this modified circuit 
into an adiabatic computation whose output is e/2 close in ^-norm to \rf) and whose running time 
is f(L',e/2). Since the trace distance between two pure states is bounded from above by the I2- 
distance (see, e.g., (Hj), we obtain that the output of the adiabatic computation is also e/2 close 
in trace distance to \rj)(rj\. In addition, it is easy to check that after we trace out the clock qubits 
from \rj), we are left with a state that is e/2 close in trace distance to the final state of the circuit. 
We complete the proof by applying the triangle inequality. ■ 

We can now apply this lemma on the result of Claim ETUI This completes the proof of Theorem 
13.11 with the running time being 0(e~^ +35 ^ L 5+2S ). 

3.1.4 Spectral Gap 

In the previous subsections, we presented a Hamiltonian H(s) and showed that inside a preserved 
subspace Sq it has a non-negligible spectral gap. This was enough for the adiabatic algorithm since 
the entire adiabatic evolution is performed inside this subspace. 

In this subsection, we show that the spectral gap of H(s) in the entire Hilbert space is also 
non-negligible. The purpose of this result is twofold. First, the existence of a non-negligible spectral 
gap in the entire Hilbert space might have some relevance when dealing with adiabatic computation 
in the presence of noise (see, e.g., [H]). Second, the techniques that we use here, and in particular 
Lemma b.121 are used again in the next subsection, but are easier to present in this simpler context. 

Lemma 3.11 For allQ<s<l, A(H(s)) = Q(L~ 3 ). 

Proof: Let S be the subspace of dimension {L + 1) ■ 2 n spanned by all legal clock states. Observe 
that S is preserved by H(s), i.e., H(s)(S) C S. Hence, the eigenstates of H(s) belong either to S 
or to its orthogonal subspace S ± . We can therefore analyze the spectrum of H$(s) and of H s ±(s) 
separately. 

First, due to the term -ff c iock and the fact that all other terms are positive semidefinite, the 
ground energy of H s ±(s) is at least 1. Second, as we will show next using Lemma l3.121 the spectral 
gap of Hs(s) is ri(L~ 3 ). To establish the same spectral gap for H(s), it is enough to show that 
the ground energy of Hg(s) is smaller than ^, which would mean that the spectral gap of H(s) is 
exactly that of Hs(s). Indeed, observe that 

( 7 o|#s00| 7 o) = (7o|#5 (s)|7o) = s/2 < 1/2 
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where the first equality holds because I70) £ So and the second follows from Equations [7| and El 
Therefore, the smallest eigenvalue of Hs(s) is bounded from above by 1/2. ■ 

Lemma 3.12 Let S denote the subspace spanned by all legal clock states. Then the ground state 
of H s (0) is I70), and that of H s (l) is \n). Moreover, for all < s < 1, A(H s (s)) = Q(L~ 3 ). 

Proof: We can write S as the direct sum of 2 n orthogonal subspaces So, Si, . . . , <S2™-1, defined as 
follows. For < j < 2 n - 1 and < I < L define |-y|> := \a j (£)) <g> \l e L - £ ), where \a j (£)) is the 
state of the quantum circuit at time £ if the input state corresponds to the binary representation j. 
Note that (7°) = \ The space Sj is spanned by {|7q), . . . , 17^)}- It is easy to check the following 
claim (see Figure EJ). 

Claim 3.13 The Hamiltonian Hs(s) is block diagonal in the Sj 's. 




H, 



Figure 2: Hs(s) is block diagonal. 

By Claims 13.21 13.31 13.51 and 13.131 it suffices to argue that the ground energy of Hs j (s) for any 

j 7^ is larger than the ground energy of Hs (s) by at least f2(l/L 3 ). Essentially, this follows from 
the penalty given by the term -ffi npu t to nonzero input states. The proof, however, is slightly subtle 
since Hi nput assigns a penalty only to states (7^) with £ = 0. 
Notice that 

HSj(s) = H So (s) + H Sj i n p U t. 

Moreover, for 1 < j < 2 n — 1, ifs^input is diagonal, with its top-left element at least 1 (it actually 
equals the number of l's in the binary representation of j) and all other diagonal elements zero. 
Hence, if we define M as 



M 



(I 


\ 



\ 


J 



then Hs^input — M is positive definite and therefore we can lower bound the ground energy of 
Hsj (s) with the ground energy of H$ (s) + M. For this, we apply the following geometrical lemma 
by Kitaev (Lemma 14.4 in [22]). 



Lemma 3.14 Let H\,H<z be two Hamiltonians with ground energies 01,02, respectively. Suppose 
that for both Hamiltonians the difference between the energy of the (possibly degenerate) ground 
space and the next highest eigenvalue is larger than A, and that the angle between the two ground 
spaces is 9. Then the ground energy of H\ + H2 is at least a± + 02 + 2Asin 2 (#/2). 



16 



We now apply this lemma to H$ (s) and M. By Claim RT31 the spectral gap of H$ (s) is Q(l/L 2 ). 
The spectral gap of M is clearly 1. Moreover, using Claim l3~Tl we obtain that the angle between 
the two ground spaces satisfies cos(6>) < 1 — 1/L by the monotonicity property of the ground state 
of Hs Q (s) (see Claim l3~T)) . It follows that the ground energy of Hs^s) is higher by at least f2(l/L 3 ) 
than that of Hg Q (s). ■ 

Remark: Notice that we only used the following properties of ^nput: its restriction to So is and 
its restriction to Sj for any j ^ is a diagonal matrix in the basis |To)> ■ • • j |Tl) whose top-left 
entry is at least 1 and all other entries are non-negative. This observation will be useful in Section 
11 



3.2 Three-local Hamiltonian 

We now show that adiabatic computation with 3-local Hamiltonians is sufficient to simulate stan- 
dard quantum computations. 

Theorem 3.15 Given a quantum circuit on n qubits with L two-qubit gates implementing a unitary 
U, and e > 0, there exists a 3-local adiabatic computation AC(n + L, 2, i?i n it, -fffinab e ) whose running 
time is poly(L, ~) and whose output state is e-close (in trace distance) to U\Q n ). Moreover, Hi n it 
and -fffmai can be computed by a polynomial time Turing machine. 

The proof of this theorem builds on techniques developed in previous subsections. The tech- 
niques developed in this section will be used in Section H 



3.2.1 The Hamiltonian 

Consider the Hamiltonian constructed in Subsection 13.1.11 Notice that all terms except Hg are 
already 3-local (some are even 2-local or 1-local). In order to obtain a 3-local Hamiltonian, we 
remove two clock qubits from the 5-local terms in Hi and leave only the Ith clock qubit. More 
precisely, for 1 < £ < L define 

H[ := 7® |100)(100|t W+ i - Ut® \l)(0\j - U\ ® |0)(l|, c + 1® |110)(110|?_ 1>W . 

For the boundary cases I = 1, L we define 

H[:=I® |00)(00|; j2 - U x ® |1)(0|5 - u{ ® |0)<1|S + I® |10)(10|f i2 

1 |C 

-1,L- 



h' l -.= i® |io)(io|i_ ljL -u e ® \i)(o\i -u{® |o)(i| c L + / ® |n>(ii|L 



Note that because of the terms |l)(0| c and |0)(l| c , these Hamiltonians no longer leave the 
subspace S invariant. To mend this, we assign a much larger energy penalty to illegal clock states. 
As we will see soon, this makes the lower part of the spectrum of our Hamiltonians behave essentially 
like in their restriction to S. Set J = e~ 2 L e and define 

^init — -^clockinit ~i~ -f^input "I" J ' -f^clock 
L 



#final : ~ 9 X/ H 'l + ^ in P ut + J ' -^clock- 



The Hamiltonian we use here is thus 

H'(s) = (l-s)H[ Jlit + sH^ al . 
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Essentially the same proof as that of Claim shows that I70) is a ground state of H' init . However, 
it turns out that |r/) is no longer a ground state of H' &nal (the proof of Claim does not apply 
since Hg is no longer positive semidefinite). However, as we shall see later, \rj) is very close to the 
ground state of H' &nal . 

3.2.2 The Spectral Gap 

Our first claim is that, when restricted to S, H' and H are identical. 
Claim 3.16 For any < s < 1, H s (s) = H' s (s). 

Proof: Let be the orthogonal projection on S. Then our goal is to show that ILsH(s)ILs = 
UsH' (s)Tls- The only difference between H(s) and H'(s) is the factor of J in -ff c i oc k) and that the 
Hi terms are replaced by H' t We note that -f^dock is zero. Hence, it suffices to show that for all 
1 <i <L, 

UgH^s = U s H' e U s . 
For this, observe that for any 1 < £ < L, 

n 5 |i)<opi 5 = |iW-')<i'- 1 o I '-c*- 1 >| c = n 5 |iio)(ioo||_ w+1 n 5 

and similarly for |0)(1|£. A similar statement holds for £ = 1,L with the right hand term modified 
appropriately. ■ 

Lemma EE] and Claim EH imply that A(H' s (s)) = fl(L~ 3 ). We now want to deduce from 
this a lower bound on A(H'(s)), without the restriction to S. For this we use the following claim. 
Essentially, it says that if J is large enough, then the lower part of the spectrum of H'(s) is similar 
to that of H' s (s). More precisely, it shows that the lowest eigenvalues, the second lowest eigenvalues, 
and the ground states of the two Hamiltonians are close. Intuitively, this holds since the energy 
penalty given to states in S ± , the orthogonal space to S, is very high and hence any eigenvector 
with low eigenvalue must be almost orthogonal to (and hence almost inside S). We note that 
a similar lemma was used in PO] in the context of quantum NP-complete problems. 

Lemma 3.17 Let H = H\ + H2 be the sum of two Hamiltonians operating on some Hilbert space 
7i = S + . The Hamiltonian H2 is such that S is a zero eigenspace and the eigenvectors in S 1 - 
have eigenvalue at least J > 2K where K = Let a and b be the lowest and the second lowest 

eigenvalues of H$ and let a' and b' be the corresponding quantities for H . Then the lowest eigenvalue 
of H satisfies a — j^k — a ' — a an< ^ secon d lowest eigenvalue of H satisfies b' > b — -j^k • V' 
moreover, b> a then the ground states of Hs,H respectively satisfy 

ma 2 > 1-- 



)-a)(J-2K)' 
Proof: First, we show that a' < a. Using fl^]^) = 0, 

(£\H\0 = + (Z\H 2 \0 = a 

and hence H must have an eigenvector of eigenvalue at most a. 
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We now show the lower bound on a'. We can write any unit vector \v) £ TC as \v) = ai\vi) + 
c^2 1 ^2) where |i>i) E S and \v 2 ) £ are two unit vectors and ai,a 2 are two non-negative reals 
satisfying a 2 + a 2 . = 1. Then we have, 

(v\H\v) > {v\H x \v) + Ja\ 

= (1 - al)(vi\Hi\vi) + 2aia2Re(vi\Hi\v 2 ) + a\(v 2 \Hi\v 2 ) + Ja\ 

> <ui|^i|ui) - ifa 2 , - 2ifa 2 - if a 2 , + Ja\ 

= {v^H^vx) + {J -2K)al-2Ka 2 

where we used a 2 = 1 — a 2 and ai < 1. Since (J —2K)o3, — 2Ka 2 is minimized for «2 = K/(J — 2K), 
we have 

(v\H\v) > { Vl \ Hl \ Vl ) - j^^. (10) 

We obtain the required lower bound by noting that (v i\Hi\vi) > a. 

Consider now the two-dimensional space C spanned by the two eigenvectors of H corresponding 
to a' and b' . For any unit vector \v) £ C we have {v\H\v) < b' . Hence, if C contains a vector \v) 
orthogonal to S, then we have b' > (v\H\v) > J — K > K > b and we are done. Otherwise, the 
projection of C on S must be a two-dimensional space. Being two-dimensional, this space must 
contain a vector orthogonal to Let \v) be a vector in C whose projection on S is orthogonal to 
If >■ Bv ((TIT)) , y > (v\H\v) >b - j^g, as required. 

Finally, let = |(f|OI 2 - Then we can write |£) = v^lf') + V^P^' 1 ) for some unit vector 
|£ /J -) orthogonal to Since |f) is an eigenvector of 

a = <f|tf|0 = /3<miO + (1-/3)^1^1^) 
> Pa' + (1-0)1/ 

K 2 \ , „,/ if 2 



a + (l-/3)(6-a) 



J — 2.K" 

Rearranging, we obtain the required bound. 



We can now bound the spectral gap of H'(s). 
Lemma 3.18 For all0<s<l, A(H'(s)) = n(L- 3 ). 

Proof: We apply Lemma 111 171 bv setting H 2 = J ■ H c lock and H\ to be the remaining terms such 
that H'(s) = H\ + H 2 . Note that Lemma 13.171 implies that the spectral gap of H'(s) is smaller 
than that of H' s (s) (which is 0(1/L 3 ) by Lemma by at most K 2 /(J — 2K). But it is easy to 
see that K = O(L), due to the fact that Hi consists of O(L) terms, each of constant norm. The 
result follows since 



This shows the desired bound on the spectral gap. Before we complete the proof, we must show 
that the final ground state is close to the history state. 

Lemma 3.19 The ground state of H'(l) is e-close to \rj). 
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Proof: Apply Lemma |3.17l as in the proof of Lemma 13.181 for the case s = 1. We obtain that 
the inner product squared between the ground state of H'(l) and \rj), is at least 1 — 5, with 
5 = {b _ a) fj_ 2K) = 0{L- l e 2 ), where we have used K = 0(L), J = e~ 2 L 6 , and b - a = Q(l/L 3 ) 
by Lemma 13.121 This implies that the ^-distance between the ground state of H'(l) and [77) is 
0{e/VL) < e. m 

We now complete the proof of Theorem 13.151 The adiabatic algorithm starts with I70) and 
evolves according to H'(s) for T = 6{e~ & L 7+3S ). Such a T satisfies the adiabatic condition (Equation 
121, using ||#g na i — ^initll = 0{L). By Theorem 12.11 the final state is e-close in ^-distance to the 
ground state of H' &nal . Lemma T3.19I implies that this state is e-close in ^-distance to \rj). Using 
the triangle inequality we note that the output of the adiabatic computation is 2e-close to \r]). The 
running time of this algorithm is 0{T ■ J ■ L) = 0(T • e~ 2 L 7 ) = 0(e"( 2+<5 ) L 14+3<5 ). 

We can now apply Lemma 13.101 to obtain a modified adiabatic computation whose output 

state after tracing out the clock qubits is e-close in trace distance to U\0 n ). The running time is 
( e -(i6+45) L i4+35) for any fixed 5 > o 

4 Two Local Hamiltonians on a Two-Dimensional Lattice 

In this section we prove Theorem ll.3l We simulate a given quantum circuit by an adiabatic evolution 
of a system of 6-dimensional quantum particles arranged on a two-dimensional grid. More precisely, 
we prove the following theorem: 

Theorem 4.1 Given a quantum circuit on n qubits with L two-qubit gates implementing a unitary 
U, and e > 0, there exists a 2-local adiabatic computation AC(poly(n, L), 6, -ff; n it, -^finai, e) such 
that Hi n n and H^ na \ involve only nearest neighbors on a 2- dimensional grid, its running time is 
poly(L, -), and its output (after performing a partial measurement on each particle) is e-close (in 
trace distance) to U\0 n ). Moreover, and H^ n& \ can be computed by a polynomial time Turing 

machine. 

As mentioned in the introduction, the main problem in proving this theorem, and more precisely, 
in moving to a two dimensional grid, is the notion of a clock. In the constructions of the previous 
section, the clock is represented by an additional register that counts the clock steps in unary 
representation. The terms Hi, which check the correct propagation in the Ith time step, interact 
between the £th qubit of the clock and the corresponding qubits on which Ug operates. If we want 
to restrict the interaction to nearest neighbors in two dimensions using this idea, then no matter 
how the clock qubits are arranged on the grid, we run into problems interacting the qubits with 
the corresponding clock qubits in a local way. The solution to this problem lies in the way we 
represent the clock. Instead of using an extra register, we embed the clock into the same particles 
that perform the computation by defining the notion of a shape of a state, to be defined later. We 
then create a sequence of legal shapes, and show how states can evolve from one legal shape to 
another. 

Although the construction of this section is more involved than the ones of the previous section, 
its analysis follows almost immediately from the analysis carried out in Theorem 13.151 To achieve 
this, we make sure that the Hamiltonians and some relevant subspaces are as similar as possible to 
those in the previous section. 
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4.1 Assumptions on the Input Circuit 

To simplify the construction of our adiabatic evolution, we first assume without loss of generality 
that the quantum circuit we wish to simulate has a particular layout of its gates. Namely, it consists 
of R rounds, where each round is composed of n nearest neighbor gates (some can be the identity 
gate), followed by n identity gates, as in Figure |21 More specifically, the first gate in each round 
is a one-qubit gate applied to the first qubit. For i = 2, . . . ,n, the ith gate is a two-qubit gate 
applied to qubits i — 1 and i. For i = n + 1, . . . , 2n the ith gate is an identity gate applied to the 
(2n + 1 — i)th qubit. These identity gates are included for convenience of notation. Any circuit can 
be transformed to such a form by introducing extra identity and swap gates. Let L = 2nR be the 
total number of gates in the circuit so obtained. Clearly, L is at most polynomially larger than the 
number of gates in the original circuit. 




Figure 3: The modified circuit with R = 3. 



4.2 The Particles of the Adiabatic Quantum System 

The adiabatic computation is performed on 6-dimensional particles, arranged on a two-dimensional 
square lattice with n rows and R + 1 columns. We number the rows from 1 (top) to n (bottom) 
and the columns from (left) to R (right). Columns number and 1 are used to simulate the first 
round of the circuit. Columns number 1 and 2 are used for the second round of computation, and 
so on. We denote the six internal states of a particle by |0}i I©)) |(D)i !©)> I©)> an d |®)- These 
six states are divided into four phases: the unborn phase |0)> the first phase |0), the second 
phase |(J[)), | (J)), and the dead phase |(x)). The two states in the first phase and the two states in 
the second phase correspond to computational degrees of freedom, namely to the |0) and |1) states 
of a qubit. We write |0) to denote an arbitrary state in the subspace spanned by |0) and |0). 
Similarly, |(J)) denotes a state in the space spanned by |(J)) and The phases are used to define 
the shape of the basis states. A shape of a basis state is simply an assignment of one of the four 
phases to each particle, ignoring the computational degrees of freedom inside the first and second 
phase. These shapes will be used instead of the clock states of the previous section. 

4.3 Geometrical Clock 

We now describe the way we represent clock using shapes. In the previous constructions, the space 
S of dimension 2 n (L + 1) was the ground space of the clock, i.e., the space spanned by legal clock 
states. Inside the clock register there were L + 1 legal clock states. Note that each such clock state 
can be described, essentially, in a geometric way by the "shape" of the clock particles: how many 
l's precede how many O's. 

We now describe the corresponding subspaces involved in our construction for the two dimen- 
sional case. For each < I < L, we have a 2 n -dimensional subspace corresponding to that clock 
state. Each of these L + 1 subspaces can be described by its shape, that is, a setting of one of the 
four phases to each particle. See Figure |U for an illustration with n = 6, R = 6. The six shapes 
shown correspond to clock states £ = 0,£ = An,£ = An + 3, £ = 5n + 2, £ = 6n, and £ = 2nR 
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respectively. Notice that each shape has exactly n particles in the first or second phase. Hence, the 
dimension of the subspace induced by each shape is 2 n . As £ goes from to L, the shape changes 
from that shown in Fig. 0Ji to that shown in Fig. 0J The locations at which the changes occur form 
a snake-like pattern winding down and up the lattice, following the layout of the gates in the input 
circuit Fig. [3 
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Figure 4: Legal clock states 
We now describe the legal shapes more formally. 

1. The shape corresponding to clock state £ = 2nr + k for < k < n has its r leftmost columns 
in the dead phase. The top k particles in the r + 1st column are in their second phase while 
the bottom n — k are in the first phase. Particles in the remaining R — r columns are all in 
the unborn phase. 

2. The shape corresponding to clock state £ = 2nr + n + k for 1 < k < n — 1 has, as before, its 
r leftmost columns in the dead phase. The r + 1st column has its n — k topmost particles 
in the second phase, and its remaining k particles in the dead phase. The r + 2nd column 
has its n — k topmost particles in the unborn phase and its remaining k particles in the first 
phase. All remaining particles are in the unborn phase. 

The subspace S is defined as the (L + l)2™-dimensional space spanned by all legal shapes. As 
in previous sections we partition S into 2 n subspaces Sj. Each subspace Sj is spanned by L + 1 
orthogonal states |tq), . . . , |t£), defined as follows. For each < £ < L and < j < 2 n — 1, the 
shape of W) corresponds to £. The state of the n active particles (i.e., those in either the first or 
the second phase), when read from top to bottom, corresponds to the state of the circuit after the 
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Figure 5: The initial state 



first I gates are applied to an initial state corresponding to the binary representation of j; i.e., it 
corresponds to the state Ui ■ Ug_\ ■ ■ ■ U\\j). More precisely, these particles are in a superposition 
obtained by mapping this state to the state of the n active particles in the following way: |0) to 
(or (Jf) for a second phase particle) and |1) to (or for a second phase particle). We often 
denote |7°), which corresponds to the all input, by \~fe). For example, I70) is shown in Figure 03 

With the risk of being somewhat redundant, let us now give an alternative description of the 
states |7g), . . . , \rf L )- This description is more helpful in understanding the Hamiltonians H" which 
we will define shortly. Consider a state (7^) for some t = 2rn. The n particles in the rth column 
are in their first phase and their computational degrees of freedom correspond to the state of the 
circuit's qubits at the beginning of the rth round. Particles to the left of this column are dead, 
those to the right of this column are unborn. The state |-)^ , 1 ) is obtained from |7^) by changing 
the topmost particle in the rth column to a second phase particle and applying the first gate in 
the rth round (a one-qubit gate) to its computational degrees of freedom. Next, the state |tL_ 2 ) 
is obtained from |7L_i) by changing the second particle from above in the rth column to a second 
phase particle and applying the second gate in the rth round (a two-qubit gate) to both this particle 
and the one on top of it. We continue in a similar fashion until we reach |7j +n ), in which the entire 
rth column is in the second phase. We refer to these steps as the downward stage. 

Next, let us describe the upward stage. The state |7£ +n+1 ) is obtained from \"fj +n ) by 'moving' 
the bottommost particle in the rth column one location to the right. More precisely, the bottommost 
particle changes to the dead phase and the one to the right of it changes to the first phase. The 
computational degrees of freedom are the same in both states. This corresponds to the fact that the 
n + 1st gate in a round of the circuit is the identity gate. 7 Continuing in a similar fashion, we see 
that the upwards stage ends in the state \^/j +n+n ) = l72(r+l)n) ma tches the above description 
of the first state in a round. 



4.4 The Hamiltonian 

The initial and final Hamiltonians are defined as 



TTll ._ TT" I If" 1 T _ TTll 

init ' clockinit ' input ' " ' clock 



1 L 

-^final := 2 X/ ^1 + ^input + J ' H< 



11 

clock ' 



= 1 

-2 



where J = e ■ L . These Hamiltonians are chosen to be as similar as possible to the corresponding 
Hamiltonians in previous sections. For example, -ff"iock has as its ground space the space of legal 



7 We could allow arbitrary one-qubit gates here instead of identity gates. This leads to a slightly more efficient 
construction but also to more cumbersome Hamiltonians. 
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clock states, S. As before, it allows us to essentially project all other Hamiltonians on 5, by 
assigning a large energy penalty to states with illegal shape. Also, the Hamiltonians H'J (once 
projected to S) check correct propagation from one step to the next. Other terms also serve similar 
roles as before. 

Let us start with the simplest terms. Define 

n 
i=l 

The indices indicate the row and column of the particle on which the Hamiltonian operates. This 
Hamiltonian checks that none of the particles in the leftmost column are in |0). Then, define 

fidocidnit = (J - ICDXCDI - ICD>«DI)i,i- 

This Hamiltonian checks that the top-left particle is in a |®) state. The remaining terms are 
described in the following subsections. 

4.4.1 The Clock Hamiltonian 

The shapes we define satisfy the following important property: there exists a two-local way to verify 
that a shape is legal. This allows us to define a two- local clock Hamiltonian, H" lock , whose ground 
space is exactly S, the (L + l)2"-dimensional space spanned by all legal shapes. 



Forbidden 


Guarantees that 


O0,O®,O® 


O is to the right of all other qubits 


0®,®®,®® 


® is to the left of all other qubits 


o®,®o 


O and ® are not horizontally adjacent 


00, 0®, 
®0, 0® 


only one of 0, ® per row 


O ® 
® ' ® ' ® 


only ® above ® 




' ® ' ® 


only below 


O ® 
® ' 


O and ® are not vertically adjacent 


® ® 
O ' 


no below ® and no ® below (x) 



Table 1: Local rules for basis state to be in S 



Claim 4.2 A shape is legal if and only if it contains none of the forbidden configurations of Tabled 

Proof: It is easy to check that any legal shape contains none of the forbidden configurations. For 
the other direction, consider any shape that contains none of these configurations. Observe that 
each row must be of the form ®*[®, ®]0*, that is, it starts with a sequence of zero or more ®, 
it then contains either ® or ®, and then ends with a sequence of zero or more O- Columns can 
be of three different forms. Read from top to bottom, it is either ®*®*, ®*®*, or O*®*. It is 
now easy to verify that such a shape must be one of the legal shapes. ■ 
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Using this claim, we can define a two-local nearest-neighbor Hamiltonian that guarantees a 
legal shape. For example, if the rule forbids a particle at location in state O to the left 

of a particle at location (i,j + 1) in state 0, then the corresponding term in the Hamiltonian is 
(10) ®)(0) Summing over all the forbidden configurations of Tableland over all 

relevant pairs of particles, we have 

-^clock := XT 
rSrules 

Note that the ground space of H" lock i s the (L + l)2 n -dimensional space S. 



4.4.2 The Propagation Hamiltonian 

The choice of legal shapes has the following important property: the shape of £ and that of I + 1 
differ in at most two locations. This means that for any t and j, the shape of |t«_ 1 ) and that of 
|7») differ in at most two locations. Moreover, if we consider the state of the n active particles in 
both states we see that these differ on at most two particles, namely, those on which the £th gate 
in the circuit acts. Crucially, and this is where we use our assumption on the form of the circuit 
(Figure EJ, the particle(s) on which the lib. gate acts are at the same location as the particle(s) 
whose phase changes. It is this structure that allows us to define the Hamiltonians H". These 
Hamiltonians act on two particles and 'simultaneously' advance the clock (by changing the shape) 
and advance the computational state (by modifying the state of the active particles). Since \^g) 
differs from |t^_i) in at most two adjacent lattice sites, this can be done using a two-body nearest 
neighbor Hamiltonian. 

The definition of H" depends on whether i is in the downward phase (i.e., is of the form 2rn + k 
for 1 < k < n) or in the upward phase (i.e., is of the form 2rn + n + k for 1 < k < n). We first 
define H'J for the upward phase. Assume i = 2rn + n + k for some < r < R, 1 < k < n and let 
i = n — k + lbe the row in which |t^_i) and \jg) differ. Then, 



H'J 



+ 



o 
® 



® 



® 



i+l,r 



i+l,r 



+ 



+ 



owo 

CDACD 



o 





o 



i,r+l 

i-l,r+l 
i,r+l 



0)(®,OI + l®,0)(0,OI) (i , r )(i, r+ i) 
(I®, OX®, 01 + 1®, 0)(0, OI) ( i, r )(i, r+ i) • 



The first line corresponds to changing the state |®,0) m t° I®)®)- The second line is similar 
for |0, O) and |®,®)- The purpose of the first two terms in each line is the same as that of 
|100)(I00| C and |110)(110| c in Hg from previous sections. 8 The difference is that here, to uniquely 
identify the current clock state, we need to consider particles on top of each other. The remaining 
terms in each line correspond to |100)(110| c and |100)(110| c in Hg 

8 There are other (equally good) ways to define these terms. For example, it is possible to define them so that 
they both act on the rth column. 
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For the case k = l,n, the definition is 

)\n,r + 



ttII 

rl 2rn+n+l 



rL 2rn+2n 



.X / \ X. 



.X / \ X 



+ 



1. r 

2. r 



1. r 

2. r 



OwO 
(DA® 

OWO 
®/\® 



n-l,r+l 
n,r+l 

n-l,r+l 
n,r+l 



+ ICDXCDI 



+ 10X01 



l,r+l 



l,r+l 



,O)(®,0l + l®,0)(®,OI)(n, r )(n,r+l) 
,OX®,0l + l®,0X®,OI)(n,r)(n, r+ l) 
O)(®,0l + l®,0)(0,OI)(l, r) (l, r+ l) 

O)(0,0l + l®,0)(©,OI)(i, r) (i, r+ i)- 

= 2nr + k and 



For the downward stage, H'J checks that a gate is applied correctly. For 
1 < k < n we define 



H'J:-. 







-U} 



+ 



+ 



0/\0 0/\0 



+ 



0/\0 
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+ 



+ 



0/\0 



k-l,r 
k,r 

k,r 

k+l,r 



The last two terms are meant, as before, to replace the terms |110)(110| c and |100)(100| c . Once 
again, to uniquely identify the current clock state, we need to consider particles on top of each 
other. The first term represents a Hamiltonian that acts on the two particles in positions (k, r) 
and (A; + 1, r). These particles span a 36-dimensional space. The matrix shown above is in fact the 
restriction of this Hamiltonian to the 8 dimensional space spanned by 



(recall that Ug acts on two qubits and is therefore a 4 x 4 matrix). Everywhere else in this 36 
dimensional subspace, this Hamiltonian acts trivially, i.e., is 0. 

For the case k = n we slightly modify the terms that identify the clock states, 
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For the case k = 1 we have 
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0/\0 



1. r 

2. r 



where the first term shows the restriction an operator acting on the particle (1, r) to the four 
dimensional space spanned by |0), |0), (recall that f7 2nr+1 is a one-qubit gate). 



4.5 Spectral Gap 

The analysis of the spectral gap follows almost immediately from that in Subsection 13.2.21 The 
main effort is in verifying that the restriction of each of our Hamiltonians to S is identical to the 
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restriction of the corresponding Hamiltonian in previous sections to S, when both are constructed 
according to the modified quantum circuit of Subsection l4.ll This, in fact, does not hold for H" n t , 
whose projection is not quite the same as that of H m ^ nt ] still, it is similar enough for the analysis 
in Subsection 13.2.21 to hold. 

Claim 4.3 H£ clockinit = H S:dockinit 

Proof: Both Hamiltonians are diagonal in the basis with eigenvalue for £ = and eigenvalue 
1 for any £ > 0. ■ 

Claim 4.4 For any 1 < £ < L, Hg t = H s/ . 

Proof: It is straightforward to verify that both Hamiltonians, when restricted to 5, are equal to 

EW^l + bi-iX-i-il - HH-i-il - hum 

j=0 

■ 

For H' } ' ut the situation is similar, although in this case the restriction to S is not exactly the 
same. Still, the resemblance is enough for the same analysis to hold: 

Claim 4.5 Both HsjnpxA an d ff^ input are diagonal in the basis |tJ). Moreover, the eigenvalue 
in both Hamiltonians corresponding to |7jjJ} for £ = is exactly the number of 1 's in the binary 
representation of j. 

Proof: Easy to verify. ■ 

The similarity between the two Hamiltonians breaks down as follows. While the eigenvalues cor- 
responding to for £ > are in Hs.input, those in H§ input might be positive (namely, for 

< I < n, the eigenvalue of is the number of l's in the last n — £ digits in the binary repre- 
sentation of j). Nevertheless, due to the remark at the end of Subsection 13. 1 .41 Lemma 13 . 1 21 holds 
here as well. We get: 

Lemma 4.6 For any < s < 1, H'^(s) has a spectral gap ofQ(L~ 3 ). Moreover, the ground state 

°/^5,final is 

The rest of the proof of Theorem 1 1.31 is essentially the same as in Subsection l3.2.2l By applying 
Lemma 13.171 we obtain that 

Lemma 4.7 For all < s < I, A(H"(s)) = ft(L~ 3 ). Moreover, the ground state of H"{\) is 
e-close to \rj). 

The proof is similar to that of Lemmas 13.181 and 13.191 This enables us to adiabatically generate 
the history state with exactly the same running time as in the three-local case (when the number 
of gates is that of the modified circuit of Subsection 14. 1|) . 

Finally, we would like to apply Lemma 13.101 as before. However, we cannot quite do this due 
to a technical issue: our Hilbert space is no longer a tensor product of computation qubits and 
clock qubits and tracing out the clock qubits is meaningless. Nevertheless, a minor modification of 
that lemma still applies. We first add, say, L/e identity gates to the end of the (modified) circuit. 
Now, the adiabatic computation produces a state close to the history state. We then measure the 
shape of the system without measuring the inner computational degrees of freedom. Due to the 
additional identity gates, with all but e probability, the outcome of the measurement is a shape £ 
for £ > L. If this is the case then the state of the system is such that the active particles are in the 
final state of the circuit, as desired. This completes the proof of the theorem. 
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